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Abstract 

Recently, a novel adaptive wave model for financial option pricing has been proposed 
in the form of adaptive nonlinear Schrodinger (NLS) equation [I], as a high-complexity 
alternative to the linear Black-Scholes-Merton model [21 13] . Its quantum-mechanical ba- 
sis has been elaborated in [4] . Both the solitary and shock- wave solutions of the nonlinear 
model, as well as its linear (periodic) quantum simplification are shown to successfully fit 
the Black-Scholes data, and define the financial Greeks. This initial wave model (called 
the Ivancevic option pricing model) has been further extended in [5], by providing the 
new NLS solutions in the form of rogue waves (one-rogon and two-rogon solutions). In 
this letter, I propose a new financial research program, with a goal to develop a general 
wave-type model for realistic option-pricing prediction and control. 
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1 Introduction 

Recently, a novel adaptive wave model for financial option pricing has been proposed in the 
form of adaptive nonlinear Schrodinger (NLS) equation pQ, as a high-complexity alternative 
to the linear Black-Scholes-Merton model [21 IS]- The new option-pricing model, representing 
a controlled Brownian motion, comes in two flavors: (i) for the case of constant volatility, it 
is denned by a single adaptive NLS equation, while for the case of stochastic volatility, it is 
defined by an adaptive Manakov system of two coupled NLS equations. 

The adaptive wave model with stock price s, volatility a and interest rate r has been 
formally defined in [H 0] as a complex- valued, focusing NLS equation, defining the time- 
dependent option-price wave function ip = ip(s,i), whose absolute square \ip(s,t)\ 2 represents 
the probability density function (PDF) for the option price in terms of the stock price and 
time. In natural quantum units, this adaptive (1+1)-NLS equation (called the Ivancevic 
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option pricing model)^ 



id t ip = ~ad ss ip-m 2 ^, (i = v /z T), (1) 

where /3 = /3(r,w) denotes the adaptive market-heat potential, so the term V(ip) = — /3\ijj\ 2 
represents the ^—dependent potential field. In the simplest nonadaptive scenario (5 is equal 
to the interest rate r, while in the adaptive case it depends on the set of adjustable synaptic 

weights {w l j} as: /3(r, w) = r Y^i=i w \ er f I 

2 Building the Wave Model for Option Pricing Evolution 

The development of the general wave option-pricing model can be summarized as follows. 

Firstly, in the case of low interest-rate r < 1, we have j3(r) -C 1, so V(ip) —> 0, and 
therefore equation ([1]) can be approximated by a quantum-like option wave packet, 'physically' 
associated with a free quantum particle of unit mass. This linear wave packet is a solution 
of the linear Schrddinger equation (with zero potential energy, Hamiltonian operator H and 
volatility a playing the role similar to the Planck constant) [1]: 

iadti/j = Hip, where H = ——d ss . (2) 

The general solution to ([2]) is usually formulated as a linear combination of de Broglie- 
type option plane-waves with the wave number k, linear momentum p = ok, wavelength 
Afc = 2-K/k, angular frequency ui^ = ak 2 /2, and oscillation period Z)% = 2ir/uik = Air/ak 2 - all 
comprising the option wave-packet: 

n 

jacket (s, t) = J2 CiAjV"-"**), (with n € N), (3) 

i=0 

where A is the amplitude of the option wave, the angle (ks — oo k t) = (ks — ^-t) represents 
the phase of the wave tpk with the phase velocity: v k = uot/k = ak/2. In addition, the group 
velocity of an option wave-packet is given by: v g = dco^/dk. It is related to the phase velocity 
Vk of a plane wave as: v g = v k — Xkdvk/d\k- Closely related is the center of the option wave- 
packet (the point of maximum amplitude), given by: s = tdut/dk. The option wave-packet 
has been used used in [1] to successfully fit the Black-Scholes call and put options data. 

Secondly, the NLS equation ([T]) has been exactly solved in PQ using the power series 
expansion method of Jacobi elliptic functions. Out of a series of solutions, the two most 
important ones (from dynamical perspective) are: 



1 Physically, the NLS equation JT]) describes a nonlinear wave (e.g., in optical fibers or Bose-Einstein con- 
densates) defined by the complex-valued wave function ip(s,t) of real space and time parameters. In the 
present context, the space-like variable s denotes the stock (asset) price. 
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1. The envelope shock-wave solution (or, 'dark soliton') H 

J[ks-±crt(2+k 2 )} 



-a 



V'shockO, t) = ±\l — tanh(s - akt) e 1 



and 



(4) 



2. The envelope solitary-wave solution (or, 'bright soliton'): 

VfeolitontM) = ± 



^ sech(s - akt) e ^s-\at(k 2 -i)\ _ 



(5) 



The adaptive NLS-PDFs of the combined shock-wave type ([!]) and soliton type ([5]) have 
been used in pQ to successfully fit the Black-Scholes call and put options data. Besides, the 
adaptive NLS-based Greeks (Delta, Rho, Vega, Theta and Gamma) have been defined, as 
partial derivatives of the shock- wave solution (J3]). 

Thirdly, two new wave-solutions of the NLS equation ([1]) have been provided in [5] , in the 
form of rogue wavesjl using the deformed Darboux transformation method developed in [8]. 



1. The one-rogon solution: 



a 



^lrogon(M) = a \/yj 



1 



4(1 + ace 2 t) 



l + 2a 2 (s-akt) 2 + a 2 a 4 t 2 

where a and k denote the scaling and gauge. 
2. The two-rogon solution: 



j[ks+<T/2(a 2 -k 2 )t\ 



a/3 > 0, 
(6) 



2rogon 



(s,t) 




1 + 



P 2 (s,t)+iQ 2 (s,t) 



R2(s,t) 

where P 2 ,Q2, R2 are certain polynomial functions of s and t. 



i[A S+ff /2(a»-fc»)t] j a/3>0i (7) 



2 A shock wave is a type of fast-propagating nonlinear disturbance that carries energy and can propagate 
through a medium (or, field). It is characterized by an abrupt, nearly discontinuous change in the character- 
istics of the medium. The energy of a shock wave dissipates relatively quickly with distance and its entropy 
increases. On the other hand, a soliton is a self-reinforcing nonlinear solitary wave packet that maintains its 
shape while it travels at constant speed. It is caused by a cancelation of nonlinear and dispersive effects in 
the medium (or, field). 

3 Rogue waves are also known as freak waves, monster waves, killer waves, giant waves and extreme waves. 
They are found in various media, including optical fibers [6]. The basic rogue wave solution was first presented 
by Peregrine [7j to describe the phenomenon known as Peregrine soliton (or Peregrine breather). 
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3 New Financial Research Program 

I propose a new financial research program as follows. 

Firstly, define the general adaptive wave model for option pricing evolution as a (linear) 
combination of the previously defined particular solutions to the adaptive NLS-equation (pQ). 
The five wave-components of this general model are: 

1. the linear wave packet Vpacket (s > t), given by ([3]); 

2. the shock-wave ip s hock(s,t), given by (H|); 

3. the soliton ^solitoa(M)> given by ©; 

4. the one-rogon V'lrogon^, t), given by ©; and 

5. the two rogon ^2rogoa(s, t) ©. 

Formally, the general adaptive wave model is defined by: 

V'general^, t) 

± 
± 
+ 
+ 

where Aj, (i = 1, 5) denote adaptive amplitudes of the five waves, while the other param- 
eters are defined in the previous section. 

Secondly, we need to find the most representative financial index or contemporary markets 
data that clearly show in their evolution both the efficient markets hypothesis [9j and adaptive 
markets hypothesis [ID]. Once we find such a representative data, we need to fit it using our 
general wave model ([8]) and the powerful Levenberg-Marquardt fitting algorithm. I remark 
here that, based on my empirical experience, the general wave model ([8]) is capable of fitting 
any financial data, provided we use appropriate number of fitting coefficients (see [U 2] for 
technical details). 

Once we have successfully fitted the most representative market data we will have a model 
that can be used for prediction of many possible outcomes of the current global financial storm. 
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c,c 



i(ks-u k t) 



(8) 



i=0 




i[fcs-±crt(2+fc 2 )] 



A2\l —q- tanh(s — akt) e 

P 

° - ' ' - i[fcs-icrt(fc 2 -l) 



Az^l — sech(s — akt) e 

P 



A 5 a 




a 

2ft 
a 
2/3 



4(1 + aa 2 t) 



l + 2a 2 (s-akt) 2 + a 2 aH 2 



1 + 



P2(s,t) + iQ 2 (s,t) 
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